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Abstract 

This paper aims at the most comprehensive and systematic construc- 
tion and tabulation of mechanical systems that admit a second invariant, 
quadratic in velocities, other than the Hamiltonian. The configuration 
space is in general a 2D Riemannian or pseudo-Riemannian manifold and 
the determination of its geometry is a part of the process of solution. 
Forces acting on the system include a part derived from a scalar potential 
and a part derived from a vector potential, associated with terms linear 
in velocities in the Lagrangian function of the system. The last cause 
time-irreversibility of the system. 

We construct 41 multi-parameter integrable systems of the type de- 
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scribed in the title mostly on Riemannian manifolds. They are mostly 
new and cover all previously known systems as special cases, correspond- 
ing to special values of the parameters. Those include all known cases of 
motion of a particle in the plane and all known cases in the dynamics of 
rigid body. In the last field we introduce a new integrable case related to 
Steklov's case of motion of a body in a liquid. Several new cases of motion 
in the plane, on the sphere and on the pseudo-sphere or in the hyperbolic 
plane are found as special cases. Prospective applications in mathematics 
and physics are also pointed out. 
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1 Introduction 



1.1 Historical 

The direct method for constructing integrable time-reversible mechanical sys- 
tems was initiated by Bertrand in his papers [1] and [2]. In the first one he 
posed and partially solved the problem of finding forces that should act on a 
particle in the plane from the knowledge of an additional integral of motion of 
a prescribed form. In the second paper [2] Bertrand formulated the problem of 
finding forces acting on a particle in the plane so that the Lagrangian 

L = ±(x 2 +f)-V(x,y) (1) 

admits a first integral of a prescribed form. He obtained the equations satisfied 
by the coefficients of an integral, polynomial or rational (linear to linear) in the 
velocities and the components of force acting on the system and succeeded to 
solve them in certain partial cases. For the case of quadratic integral, Bertrand 
obtained the necessary and sufficient condition is that V satisfies the partial 
differential equation 

{2axy + by + b'x - c±){V xx - V yv ) + 2[a(y 2 -x 2 ) + by-b'x + c]V xy 
+3(2ay + b)V x - 3{2ax + b')V y 
= (2) 

where a,b, ... are constants. Darboux [3] gave the generic solution of this equa- 
tion for the potential in the form 



v _t^m (3) 

or — (i 



where /, tp are arbitrary functions and a, (3 are the parameters of confocal ellipses 
and hyperbolas in the plane. The cases of degeneracy of elliptic coordinates to 
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parabolic, polar and Cartesian coordinates were considered in complete form 
only recently [4]. The last cases were also rediscovered more than once in later 
times [5] , [6] [7] . The same result was drawn from different points of view in [8] , 
where some complex cases were pointed out (see also [9] for complete description 
of possible real and complex cases) and [25] , which applied differential geometric 
concepts to the iso-energetic Jacobi metric. It turns out that namely in those 
four cases the Hamilton- Jacobi equation of the system (1) is separable [11], [12] 
(see also [13]). 

In the present paper we deal with a significant generalization of the system 
of (1), described by the Lagrangian 



in which a{j,ai,V depend only on the coordinates qi,q2- This system differs 
from (1) in two important aspects: 

1. It accommodates Riemannian two-dimensional manifolds as possible con- 
figuration spaces 1 . This allows application to various problems in dynam- 
ics involving motion of a rigid body about a fixed point and motion of a 
particle on a fixed smooth curved surface. 

2. The terms linear in velocities account for gyroscopic forces. Those are 
forces that do not do work during motion of the system. They appear 
naturally in gyroscopic systems as a result of ignoring cyclic coordinates 
by the Routh procedure. Linear terms may arise also when the system 
has some charged components moving in a stationary magnetic field. In 

x In classical mechanics we are usually interested in systems whose kinetic energy (equiv- 
alently, metric) is positive definite. Certain results below involve indefinite (pscudo- 
Ricmannian) metrics. These are included for the sake of completeness of enumeration of 
all possible cases. 
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this context the vector (ai,a 2 ) may be called "vector potential" while V 
is the scalar potential. The system (4) is usually called a system with 
velocity-dependent potential. For other circumstances see [16]. 

Birkhoff [29] raised and completely solved the problem of constructing sys- 
tems of the type (4), which admit an integral linear in velocities in addition 
to the Jacobi integral (the Hamiltonian of the system). The first cases of irre- 
versible systems in the Euclidean plane, having an integral quadratic in velocities 
were obtained by Vandervoort [14] and by Dorizzi et al [15]. 

The first systematic study of the full irreversible system of type (4) for 
existence of a quadratic integral was done in our work [16], where possible 
integrable systems were classified according to whether the normal form of the 
metric has the Euclidean form, the form of a metric on a surface of revolution or 
a generic Liouville form. Several new many-parameter families of systems were 
found. Most of the known systems of the same type were either recovered or 
shown to be special cases of the new systems, corresponding to certain particular 
values of the parameters. Some main systems were written down explicitly, but 
many others were only briefly pointed out. Some cases of integrable motions in 
the plane generalizing the special plane versions of [16] were presented in [17]. 
Cases of irreversible motion on a sphere corresponding to integrable problems 
due to Clebsch and Lyapunov in rigid body dynamics were obtained in [18]. 
A case of irreversible motion of a particle on a smooth ellipsoid corresponding 
to Clebch's case of asymmetric body in a liquid was also obtained in [19]. In 
the last two works the Lagrangian for each motion was expressed in redundant 
Cartesian coordinates on the surface. They can be simply reduced to the form 
(4) by using sphero-conic coordinates on the sphere and elliptic coordinates on 
the ellipsoid. 

After the publication of [16] few works appeared dealing with the same prob- 
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lem, but mainly investigating the case of Euclidean plane. A nearly complete 
list is composed of the following papers: [20], [21] (devoted to motion of a par- 
ticle in a rotating plane) and [23], [25], [26] dealing with particles moving in 
the plane under the action of potential and "forces with vector potential" or 
"velocity-dependent" forces. Only two works considered irreversible motions on 
curved manifolds: [24] from the point of view of separation of variables and [22] 
in which a classification of possible cases on Riemannian manifolds is tried, but 
most of the results are very special cases of those presented in [16]. 

The aim of the present paper is to make the most exhaustive enumeration 
of irreversible mechanical systems, which have an invariant quadratic in veloc- 
ities and functionally independent of the Hamiltonian of the system. We have 
constructed a total of 41 different several-parameter families of those systems. 
Some of those families or their special cases or degenerations were listed in our 
previous works but the majority are listed for the first time. We give briefly the 
very essential moments in the different methods used in the derivation of the 
families and then provide for every family the minimal information enough for 
its characterization. Usually, we give the Lagrangian an the second invariant 
and occasionally the Gaussian curvature of the configuration space, which we 
use for interpretation of suitable results as cases of motion on a plane or a space 
of constant curvature. The last section is devoted to application in the two 
model problems of motion of a rigid body under various circumstances and of 
motion of a particle under potential forces on a smooth fixed ellipsoid. 

1.2 Formulation of the problem 

It is well known that the system (4) can always be referred to some isometric 
coordinates x,y (say), in which the Lagrangian takes the form 

L = ^A (x 2 + y 2 ) + hx + l 2 y-V (5) 
6 



containing one coefficient in the kinetic energy term, instead of three in (4) . The 
pair of coefficients li,l 2 are not unique for a given system, since one can add 
the time derivative of an arbitrary function of the position to the Lagrangian. 
The Lagrangian -more precisely, the Lagrangian equations of motion- is thus 
characterized by the triplet of functions (A , CIq, V), where f2 = ^ — 
The system (5) transforms under conformal mapping 

x + iy = z(( = £ + irj) (6) 

to 

L=\k{(? +r t 2 ) + U + l' 2 i 1 -V (7) 

where A =| | 2 A ,O=| | 2 fio- The last system admits the Jacobi integral 
of motion (the Hamiltonian expressed in the state variables) 
1 • 2 

H = -A(£ + r/ 2 ) + V = K (arbitrary constant.) (8) 

Suppose that the system (5) admits an integral of general quadratic form 

/ = i(6ni 2 + 2b 12 xy + b 22 f) + b lX + b 2 y + b (9) 

In our previous work [16] we have shown that for a system of the type under 
consideration one can always find a conformal mapping of the plane and a time 
transformation of the type 

dt = Adr (10) 

such that the Lagrangian, Jacobi's integral and the complementary quadratic 
integral take, with the proper use of (10), the form [16]: 



L = l^+^ + Pf-Q^ + U (11) 
H = l -(e+V a )-U (12) 

/ = 1^ + p£ + q v i + r (13) 



where primes denote differentiation with respect to r, 

P = <j>„,Q = 4>t, (14a) 

U = A(K-V), (14b) 

R = f'nPdrj- (SlQ + U ( ) v = Vo dS, (14c) 
•'■no 

& = hi + trm ( 14d ) 
and A, <j>, V satisfy the over-determined system of four equations 

A 5v = 0, (15a) 

4> v At + ^A„ + 2A<% ) = 0, (15b) 

+ W = °- ( 15c ) 



[^ + ^](W-(A^ = (ir„l) 



The equations of motion corresponding to the Lagrangian (11) acquire the 
simplest form 

In the original (natural) time variable the Lagrangian and integrals take the 
form 

L = ^A(f + i 1 2 ) + Pi-Qr 1 + K -V (17a) 
H = ^A(f + i] 2 ) + V-K = (17b) 
I = \^ 2 k 2 + A(PC + Qfi) + R (17c) 

In the last form of the integral (17c) the function R may involve the energy- 
parameter (more precisely, Jacobi's parameter) if as a linear multiplier in cer- 
tain terms. This is obvious from the way of construction. This parameter is 
interpreted as its numerical value on a given motion, or may be substituted by 
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its expression as a function of the state variables from (17b), resulting in an 
unconditional integral in the state space. 

The function A, the solution of equation (15a) completely characterizes the 
metric on the configuration space of the system (17a). This gives a natural basis 
for a classification of all systems of the type under consideration, depending on 
whether A has one of three possible forms: 

1) The factor A is a constant and, without loss of generality, one can set 

A = 1 (18) 

In that case the configuration space has Euclidean metric. It may be interpreted 
as a plane or a developable surface. 

2) The factor A depends only on one variable 

A = n( V ) (19) 

In that case the metric on the configuration space has the structure of a metric 
on a surface of revolution. This does not mean that one can always realize 
this metric on a surface of revolution, but this may be done in certain special 
circumstances. 

3) A has the generic form of a Liouville metric 

A = A(0 - /ifa) (20) 

The solution of equations (15b-15d) will depend mainly on the form of chosen 
solution for A. The three cases will be dealt with separately in the following 
sections. 

2 Time-reversible plane systems: The Bertrand- 
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Darboux project 

When the function = equations (16) become time-reversible and one can 
take <p = 0, so that in virtue of (15a, 15d) the Lagrangian can be written as 

This is the general two-dimensional system of Liouville type, and it can be put 
in the common form. 

It involves four arbitrary functions of each of one variable each. Separation of 
variables in this case is obvious. 

Many cases of physical interest can be identified in the form of the Liouville 
system (21). 

If the original configuration space of (5) is the ordinary Euclidean plane, then 
A = 1 and we obtain a different and very compact formulation of the Bertrand- 
Darboux problem presented in the introduction. Solutions of this problem are 
systems of the form (21) with the function z(£) in the transformation (6) satis- 
fying the equation 



A 7 

A=|^| 2 =A(0-M*7) 
It was shown in [16] that this equation has the general solution 

C=[- dZ (23) 

in which a is a non-negative constant and the arbitrary constants/3, 7 can be 
made real by a suitable choice of axes in the z— plane. The generic case of two 
distinct roots defines the elliptic coordinates in the z— plane, reproduces Dar- 
boux's result. The degenerate cases of two equal roots or with one or both roots 
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coalescing with the point at infinity define the polar, parabolic and Cartesian 
coordinates in the plane. This gives a much simpler and transparent way to the 
completion of the Bertrand-Darboux project than that relying on the cases of 
solution of Bertrand's equation (2). 



3 The first type of irreversible systems - The 
metric is an Euclidean plane one (System 1) 

We begin the classification by considering the first and simplest case, when 
A =const.= 1. It is quite straightforward to solve the system (15) as in [16] and 
construct the Lagrangian in the form 

L = l -(x 2 +y 2 ) + ^T)x-^F{X)v 

+ \a(X - fif + + b 2 ){X - fi) 2 +K 1 (\-n) + h (24) 



where 



J jm ,v ~l 



F(X) = aX 3 + hX 2 +aX + di, 

G(fi) = -a^i 3 + b 2 n 2 + c 2 n + d 2 (25) 

and a, bi,b 2 , C\,C2, di,d 2 ,Ki,h are arbitrary constants. This is equivalent to the 
less explicit result provided in [16]. The equations of motion have the form 

dV dV 
x + fly = - — ,y-Clx = -— (26) 



where 



= i [3a(A 2 - M 2 ) + 2(6! A + b 2i i) + a + c 2 (27) 
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They admit the complementary quadratic integral 

/ = \x 2 + ^G[iT)x + VF(A)y 

~[a{\ - fi) 2 (2X + ft) + (ci + c 2 )(A - /x) + b 2 (\ 2 - A* 2 )] - M(A + 

The system with the Lagrangian (24) describes the motion of a particle in the 
plane, with arbitrary Jacobi's constant h. An equivalent of this system was 
pointed out in [15] and discussed also in [16]. The Lagrangian can be explicitly 
expressed in terms of the Cartesian coordinates x, y using Weierstrass' elliptic 
functions of two independent sets of invariants. Different possible real interpre- 
tations are discussed in [16] in terms of Jacobi's elliptic functions. Except some 
degenerations, the potential V and the gyroscopic functions f2 are periodic in 
both x, y directions. The generic system (24) may be used as a model for the 
study of motion of an electrically charged particle in the plane under the action 
of an electric field in the plane and a magnetic field orthogonal to it. Both fields 
may be regarded as resulting from a spatially periodic distribution of sources. 

The simplest periodic case, when elliptic functions degenerate into trigono- 
metric ones is characterized by 

V = i(a 2 + /3 2 )[Acos(ax) + B cos(f3y)} 2 - k[Acos(ax) + Bcos(f3y)] 

= a 2 Acos(ax) + /3 2 B cos(/3y) 

1 = x 2 - y 2 + 4:fiBsm(l3y)x + 4:aAsin(ax)y - 2k[Acos(ax) + B cos(/3y)j 

+ L4cos(aa;) - B cos(/3y)][A(3a 2 + j3 2 ) cos(ax) + B(a 2 + 3/3 2 ) cos(fii^) 

Even this case is a quite complicated integrable one. The following two graphs 
1, 2 show the functions V and SI on an area in x 47r, containing four primitive 
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periodic cells of the plane in the special case a = fi = A = B = k = l. 



Figure 1 :V(x, y) the electric potential in the plane of motion. 



Figure 2: Q(x, y) the magnetic field orthogonal to the plane. 
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4 Second type 

The problem of constructing systems of this type was considered in [16], where 
only few cases were given in detail. Here we give a more systematic treatment 
and the final results cover a much wider hierarchy of possible different cases, 
some of which were overlooked in [16]. 

In fact, putting A = /j(r?) in (15) we obtain the general solution of (15b) 
as tp — + N (77), where M, N are arbitrary functions. Also, from (15c) we 
deduce that V must have the structure 

V = V^),i> = ^JlM(0+v(r ] ) (30) 

where v is related to N by the condition u'(r?) = —fiN'(r]). In order to make 
the remaining equation (15d) more tractable we have used A = M(£) and \x as 
variables instead of £, 77, and introduced the notation 

( ^ )2 = F(A) ' ( f f=G{ji) (31) 
where F and G are functions to be determined. Now we can write 

r-jx-f^ (32 ' 

and equation (15d) becomes a consistence condition which has to be satisfied by 
four functions in one variable each V(if>), F(X), G(p) and v(X). The last equation 
can be shown to have the structure 

/^"(V) + \v'W) = S (33) 

To solve this condition, i.e. to find compatible choices of the four functions we 
have applied two ways that led to different findings. 

4.1 The first way: (33) is a differential equation 
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For (33) to be a differential equation in V(^>), the coefficient fitjj^ and the 
function S must be functions of the single variable ij). The first implies that 
v(/j) = k^/Jl. Turning to the function S we have found that k = and also 

F{\) = A 2 (a 3 A 3 + a 2 A 2 +aiA + a ) 

G(n) = n 2 (b 3 n 3 + b 2 fi 2 + b lf i-a ) (34) 

Returning to equation (33) we find 



l lu .1.2 , fl 3^ , A [JI 



^) = --(6 3 V 2 + 4) + -2, V>=a/t (35) 



By this information we are now able to construct the system we are looking 
for. We call it System 2. It has the Lagrangian 



2 .-.2 



1 A , /i 



2^ L F(A) G(m) j 



+ 1 J / F(X)G(fi) r AA MM 



2 V A 3 m 3 ^(M) 



^4A t2 1 , A 2 m 
M 8 m A 



+ J 2 o(a 3 ^+&3T) (36) 



4.2 The second way: (33) is a consistency condition: 

We deal with (33) assuming for V(if>) a preassigned form like a polynomial and 
try to find compatible combinations of coefficients that help separate equations 
for F, G and v. The result found in this way can be formulated as in the following 



Theorem 1 
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The system with the Lagrangian 

• 2 

1 /iA jj, 



L 



2 a + ai\ - c A 2 4^(c 3 it 3 + c 2y u 2 + ci/x + c ) 

-[— + 2v'(fi)} ^ C3 ^ 3 + C2 ^ + ClM + C ° A 
v 7 ^ V a o + fl iA — c A 2 

Jy ao + aiA — coA 2 /i 

2/1 3 / 2 C 3 /X 3 + C 2 /X 2 + Ci/X + c 

- [JV7^A + w(/i)] + ^c 3 [JVmA + v(p)} 2 (37) 

where Cj,di,b,J are arbitrary constants and is the solution of the linear 
differential equation 

8(c 3 it 3 + c 2 ^ 2 + ci/i + co)v"'((j,) + 12(3c 3 /i 2 + 2c 2 /x + ci)u"(/i) 

+6(3c 3M + C2)«'(m) - 3c 3 w( A1 ) = -36 - |^ (38) 
admits the second quadratic integral 



1 fj, 2 X 



■ 2 



2 



2 ao + aiA — coA 

r ^A c 3 /i 3 + c 2 /i 2 + ci/x +~cq" ■ 

-M[-= + 2w (/i)] = A 

Voo + a iA — CoA 

J/i\/ a + aiA — c A 2 /i 



2/U 3 / 2 sjc-ip? + c 2 fi 2 + ci/x + 



^3/2 



d/i 



+ i J 2 A[Ait(c 3 ^ + c 2 ) - ^] - i Jai y 
+ JAVM{4(c 3 /i 3 + c 2 /i 2 + ci/x + co)u"(/i) 
+2(3c 3 it 2 + 2c 2 /i + ci)w'(/i) - C3(j,v(fi) + b/i} 
+2w'(ii) 2 (c 3 ii 3 + c 2 Ai 2 +ci/i + co) (39) 

When c 7^ the above theorem gives the same result as in [16], as it is 
possible in this case to eliminate a\ by a suitable shift in the variable A. But 
when Co = the presence of a\ is significant and leads to seven new cases in 
dependence on the combination of roots of G(/i). The complete hierarchy is 
given in the next subsection. 
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4.2.1 Classification 

For the complete determination of the Lagrangian (37) and the integral (39) we 
need only the solutions of (38) for v(p). This turned out to be depending on the 
distribution of roots of the polynomial G(fj,)/^i 2 = c 3 ^ 3 + c 2 ^i 2 + Cifi + c , which 
we denote by /x 1 ,^ 2 ,^ 3 . We now proceed to enumerate the fourteen possible 
cases. For each case we give its conditions and the corresponding solutionw(/i). 
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Conditions on coefficients Roots 


v{n) 


3 


c 3 ^0 


c + 


M 3 7^ M2 + Ml 7^ 


+K 2 ^//j, - M2 + -&3\/M - M 3 


4 


M 3 = M 2 7^ Mi 7^ 


^ + K ly //j, - Mi + ^2VM - M2 
+^ 3 /Vm - \x 2 


5 


M 3 = M2 = Mi 7^ 


£ + K ly //j, - Mi + ^2/VM - Mi 
+^ 3 /(M - Mi) 3/2 


6 


c = 0, ci 7^ 


M 3 7^ M2 7^ Mi = 


a' +4 C3 /;; 3 ^+^2Vm m 2 

+K 3 yjn - M3 


7 


c = 0, ci 7^ 


M 3 = M2 7^ Mi = 


£+4S37^ +A 2VM-M2 

+ *3 


8 


co = ci = 0, c 2 7^ 


M 3 7^ M 2 = Mi = 


i + J* *sv7* M 3 


9 


c = ci = c 2 = 


M 3 = M2 = Mi = 


b 1 Ki , K2 1 Jai 
c 3 " 1 " v^ 7 M 3/2 3 2c 3 M 5/2 


10 


c 3 = 0, c 2 ^ 


co ^0 


M 2 7^ Mi 7^ 


-|£ + ^iVM - Mi + ^2VM - M2 


11 


c = 


M 2 = Mi 7^ 


-277 + A iVM-Mi+^2/VM-Mi 


12 


M 2 + Mi = 


-2^-I^+WM-M2 


13 


co = ci = 


M 2 = Mi = 


6^ | ifi , Jai 
2c 2 + + 16c 2 ^ 3 / 2 


14 


c 3 = c 2 = 
ci ^0 


co ^0 






15 


co = 






16 


C3 = C 2 

= ci = 


c 7^0 




-£+^M 2 + ^iM 



Table I 



The Gaussian curvature of the configuration space of (37) is 

k = -(3c 3 ^ 2 + 2c 2 fi + ci) 

= -^'(M) (40) 



The variable [i should be restricted to an interval [a, b] (say), where a, b e 
{— 00,0,00,^,^2,^3} on which -F(m) > 0- The Gaussian curvature k on the 
configuration manifold is of variable sign in general. But it is possible to have 
manifolds of curvature of definite sign. An example is when one of the two 
bounds is taken to be the point at infinity. We note also that in order that 
the Lagrangian was real, the coefficients Ki should be taken real, imaginary or 
including a complex conjugate pair, depending on the nature and relative order 
of the roots of F. From (40) we also conclude that only case 16 lives on a flat 
manifold and only cases 14, 15 has a manifold of constant non-zero curvature. 
In the last two cases the system describes either: 

1. a motion on a sphere (when ci < 0), which can be recognized as the 
axisymmetric version of Clebsch's integrable case of motion in a liquid of 
a body with a spherical inertia tensor. The full Clebsch case is obtained 
as a special version of another new integrable system in 6.2.1. 

2. or a motion on a pseudo-sphere or in the hyperbolic plane H 2 when c\ > 0. 

4.2.2 Example: A case of motion on a Riemannian manifold of pos- 
itive Gaussian curvature 

In order that the Lagrangian (37) represents a real system, the kinetic energy 
should be positive definite. In case 6 this implies that the variable /U should be 
confined to an interval of the non-negative half of the /i-axis on which c 2 [i 2 +cifi+ 
co > 0. Various possible choices of the interval lead to different configuration 
manifolds. We give an example of motion on a convex manifold: 

In case 10 of table I, let c (= c 2 fi 1 fi 2 ) = 1, ai = 0, a = a 2 ,0 < 11 < /Uj < fi 2 , 
The substitution A = a cos tp,fi — [i 1 sn 2 u, K x = iki,K 2 — ik 2 (^ — [i 1 = 
iy/JT[cnu, ^J[i — [i 2 = iy/JI^dnu), where the modulus of the elliptic functions 
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k = . ^ transforms the Lagrangian to the form 

V ^2 



1 ..o 2 -2\ a sin os 
-^{u + sn u G9 + ^— ? 



. — r , , . , , ,b acosip . 

+V c 2[« ; i\/^2 emu + /^y^cn it + vMi7*2( F= ) cnwdnwjx 

''.! x //jsiw/ 

-^a^snu cos 93 — sn 2 u — kiy/p^cnu — ki^fp^ dnu] (41) 

2C2 

The Gaussian curvature of the metric of this system is strictly positive 



k = c 2 ( y u 1 cn 2 u + p 2 u ) > £2(^2 — Mi) > (42) 

It is possible to realize the Riemannian metric of the last system on a closed 
surface of revolution. In fact, if (p(u),ip, z (u)) is a current point of the surface 
in cylindrical coordinates, then we have 

dp 2 + dz 2 + p 2 d<p 2 = p^du 2 + sn 2 u dip 2 ) (43) 

Comparing both sides we find that 

P = VMi snu (44) 
and hence we get an equation for z as 

2 i 2 / \ 2 
p l cn iidn u +(^ — ' = Mi 

which, on separation and integration, gives 

z = VJH / \/l + & 2 — fc 2 sn 2 u sn udu (45) 
Jo 

Noting that this is a periodic function with period 4K(k), we deduce that the 
surface thus constructed is a closed one. This family of surfaces of revolution 
depends on two parameters p l and k. The first can be absorbed by a suitable 
scaling, but k changes the form of the surface. Figure 3 shows this surface for 
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Figure 3: The surface of revolution described by (44) and(45) for k = 

5 The third type - Systems with a Liouville-type 
metric 

5.1 General considerations 

Now we consider the solution of the system of equations (15) in the generic case 
when A'(£)^'(?7) ^ 0. It will be easier in this case to use A,/z as variables. Let 



A' 2 (0 = F(A), M /2 (r?)=G( M ) 



(46) 



so that 




(47) 
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Equations (15b-15d) can be written as 



(px-cp^- 2(A - M )0a m = (48a) 
+ ^ = (48b) 

-y/F(\)G(ri[(\-»)V] Xll = (48c) 

One of those equations can be readily integrated. In fact, the equation of the 
characteristics of (48b) for V is 

^ x d\ - 0„d/i = (49) 

Multiplying this equation by (X — fi) one can verify that the resulting expression 
is a total differential, in virtue of (48a), so that one can write 

# = (A-/i)(0 A dA-^d/i) (50) 

Thus, the potential V can be expressed as V — V(ip), where ip is related to 4> 
by the relations 

i>\ = ( A - m)</>a, ^ = -(a-m)^ (51) 

Since only the derivatives of <ft enter in (48a) and (48c), one can completely 
eliminate </> from those equations. This gives for if) the equation 

^ x - ^ + 2(A - M )VV = (52) 

-(A - M ) v /F(A)G(m)^"(V)V'a^ + 3^'W#a„] 
= (53) 
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The Lagrangian describing the system (in the fictitious time) may now be 
written as 

l r A' 2 m' 2 , 



L = 



2 l F(\) G(/i) J 

1 V^)G(M)[4^' + 7^Ml 



(A-/i) v w V ^ L F(A) G(/x)' 
+(A-/i)[A--V(^)] (54) 



and in the natural time it takes the form 

' 2 2 
1 \r A /i 



2^ \F(A) G( P ) S 

1 ■V^WgM[4ttA + -^a] 



(A — m) w v " ;l F(A) G(/x)' 

(55) 

In certain circumstances we need also the expression 

° = F(A) 9A ( A^ ) - G[jl W —J + 2(A—^ (56) 

The process of constructing time-irreversible mechanical systems with a 
quadratic integral is now reduced to the simultaneous solution of the pair of 
equations (52) and (53), which involve the function tp(X, it) together with three 
functions of one variable each: F(X), G(p) andV '(ip) . It is to be noted also that 
the three one-variable functions occur only linearly in the last equation. This 
becomes obvious if we write down this equation in the expanded form 



+2(A - n)\y" + 3v'(V)Vv] 



= (57) 
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5.2 Solving the equation for ip 

Equation (52) and (48a) were given an interesting interpretation in [16]. We 
provide it here and use it in a more systematic way to construct simple solutions 
of that equation, future candidates of satisfying (53). In fact, if we introduce 
new variables p = A — p and z = i(X + p) then (48a), (52) and (51) reduce to 

<t> PP + -4> P + <f>zz = ( 58a ) 
VW-j^p + V^ = (58b) 
HV)* = -p4> p , (-iip) P = P<t> z (58c) 

Imagine z as the axis of cylindrical coordinates in a three-dimensional complex 
space and p is the radial distance of the current point from that axis. Equation 
(58a) is just Laplace's equation satisfied by 4>— the velocity potential of a virtual 
flow of an ideal incompressible fluid, symmetric around the z axis, (58b) is the 
equation for tp— Stokes stream function of that flow and (58c) are the well-known 
relations between the two functions (see e.g. [41]). In [16] we have utilized some 
known axisymmctric hydrodynamic flows to construct certain solutions of (52). 

Here we provide a systematic way to construct solutions of (52) and (48a). 
It is well known that (58a) admits the spherically symmetric solution <j) = j. = 
, 1 representing a source at the origin. For a source on the z— axis at a 
distance zq from the origin the velocity potential becomes 



1 



y/zl- 2z z + z 2 + p 2 
1 



\J z 2 - 2iz a (X + p) - 4\p 
This function has two types of expansions: 



(59) 
(60) 



= E J*T*n (61) 



2=0 Z 
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for large z and 

00 eh 

'■S^ (62 » 

for small z a , where <f> n is a homogeneous polynomial of degree n in A, fi and 
{k n },{k' n } are numerical (real or imaginary) coefficients, which can be used to 
give the polynomials <j> n the most convenient form. Since the expansion is valid 
for arbitrary (large or small) values of zo, we conclude that the two infinite 
sequences of functions {<fi n }, { } are solutions of equation (48a). It is 
noteworthy to mention that those functions correspond to the familiar sequences 
of axisymmetric solutions r™P„(cos 9), r - (" +1 )p„(cos 9) of Laplace's equation in 
spherical coordinates. 

We now follow a parallel line for solving equation (52). The spherically sym- 



metric solution of Stokes equation (58b) is easily seen to be ip = r = \J z 1 + p 2 . 
Displacing the source as indicated above we get Stokes function 

i\) = y / z2-2iz (A + /i) -4A^ (63) 

For this function we obtain the two expansions 

* = E ^» (64) 
n=0 Z 

and 

00 , 

V,= V^E^0 7^ (65) 
thus giving the two infinite sequences of solutions of Stokes equation 58b: 
{ip n ,n = 0, 00}, {-</>„ = -n^fer , n = l,...,oo}. The first few non-constant 
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ip n 's are 



A = 1 

ipi = A + /i 

V<2 - (A-M) 2 

1>3 = (A - i*) 2 (A + /i) 

V> 4 = (A-/i) 2 (5A 2 + 5m 2 + 6A/i) 

V> 5 = (A-/i) 2 (A + A1 )(7A 2 + 7/i 2 + 2A A1 ) (66) 



A + it (A -it) 2 (A- i u) 2 (A + ) u) (A- i u) 2 (5A 2 + 5^ 2 + 6AAi) 



Finally, we note that equation (52) is invariant with respect to shifting (A, /j.) 
to (A — v,fj, — v) for all constant u, so that the sequences of solutions (66,67) 
remain valid after this shift. We shall use this freedom in certain circumstances. 

5.3 Adapting a solution 

A possible way of solving (52-57) is choosing a solution /x) for the linear 
equation (52) and inserting this solution in (57), which will serve as a consistency 
condition. A working choice ip(\,fi) should allow the existence of F,G,V, so 
that this equation is satisfied. It is natural that some solutions of (52) will be 
accepted in (57), while others will be excluded. Till now, we have not been able 
to separate any part of the last equation to operate independently on one or 
two functions. It remains a matter of guess and trial to find possible cases. It 
is not possible in this way to claim that this way will lead to construction of all 
possible cases of the type under consideration. However, it has lead to several 
rich families of systems in one of the ways described below. 



The first few of the second sequence {V>„} are 





(A/x)V2 



(67) 
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5.3.1 The first procedure: 

In general, equation (57), which may be written in the form 

v"^) + ^-3V'(iP) = N(\,ri (68) 



where 

N = { [F"(A) + G»] 



(A-m) 2 



+2{[ (X^F ]aaF(a) _ [(T^WG(m)}}/[2(a - MMW) 

is a consistency condition which should be satisfied by the four functions F(X), G(n), V(ip), 
and ip(\,[i), of which only ip is subject to another equation, namely, (52). An 
important situation is met when (68) becomes a differential equation in V(tp). 
For this to happen, the coefficient of V' as well as the function N should be 
functions of tp. That is 

A^-AVAa = (70b) 

Thus, for the present procedure to work, ip must satisfy simultaneously the 
linear second-order equation (52) and the nonlinear third-order equation (70a). 

We now begin with trying a solution of (52) in the form of one of the functions 
ip n and %j) n from (66,67). It turned out that only four functions ip x , ip 2 , Vh> 
satisfy (70a). Moreover, we have tried solutions in the form of linear combina- 
tions 

5 
i=l 

This yielded the above four functions and two more functions ip l ± = 
(\/X± v 7 ^) 2 ' so that the present procedure can be applied for the six cases: 
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J(A + ft) 
J(\-rf 

jh+±L 

j(Va + vm) 2 
j(Va- vm) 2 

Table II 

Those cases will be discussed in the next section. For each case: 

a) The function ip is substituted in (70b) and separation of variables is 
achieved after affecting a differential operator. 

b) F(\), G(p) are determined. It turned out that F is rational in case 17 and 
polynomial in all other cases and also G(p) = —F (p) in all cases except 
19 and 20. 

c) F(X), G(p) are substituted in (69) and the expression N is expressed as a 
function of ip only. 

d) N(i[)) is inserted in (68), and the last is solved for V{ip). 

e) The final form of V is obtained by expressing ip again in terms of A, \i. 

Some of the indicated steps are quite cumbersome, even with the use of com- 
puter algebra packages. For space saving most of those details are not reflected 
here. We shall give only the data necessary to unambiguously reproduce every 
case. 



case 

17 

18 

19 

20 

21 

22 



28 



5.3.2 The second procedure: 

The last procedure determined cases in which the potential V is a solution of a 
second order ODE (68) depending linearly on one significant arbitrary integra- 
tion constant. The other being an immaterial additive constant. This does not 
exclude the existence of other potentials not containing this arbitrary constant. 
It may of course depend on the parameters appearing in other functions ip, F, G. 
For such cases is not required to satisfy (70a-70b). One may try expressing 
V(ip) in the form of a polynomial (say) 

JV 

V = J2^\ (72) 

i=i 

choose a solution tp of (52) and insert this expression in (57) and try to find com- 
binations of parameters that lead to separation of the two unknown functions 
F(A),G(/x). Obviously the outcome of this method will depend on the ability 
to perform necessary operations on functions V, ip involving as much coefficients 
as possible. 

We have made three trials: 

N 

A) For tp = ^ Jitp i the separation resulted in 12th-degree polynomials for 
F and G, but the final analysis led to 6 cases 

(a) Two special versions of cases 18 and 19 of table II. 

(b) ip is a combination of ip 1 ,ip 2 ,V2 = 0,F,G are cubic and G([i) = 
-Ffr). 

(c) ip is a combination of tp 1 , ip 3 , V2 = 0, F, G are quadratic and G(fi) = 
-F(»). 

(d) ip is a combination of ip 1 , ...,ip 4 ,V2 — 0,F,G are linear and G(p) = 
-F(p). 
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(e) V is a combination of ip 1 , ...,ip 5 , v 2 = 0, F(fi) = 1 = — G(/z). 

It was noted that the presence of any of the few next with ii > 5 
led to the inconsistent result G(/z) = F(/i) = 0. 

5 

B) For ip — Jj^ we have been lead to 5 cases: 

i=l 

(a) A special version of case 17 of table II. 

(b) V is a combination of r/> 2) = 0, F, G are of 5th degree with two 
equal roots and G(/x) = —F(fj,). 

(c) V is a combination of ...,^> 3 ,v\ = 0, F, G are of 5th degree with 
three equal roots and G(/z) = — F(/z). 

(d) ^ is a combination of ...,^> 4 ,v\ = 0,F,G are of 5th degree with 
four equal roots and G(p) = —F(fi). 

(e) V is a combination of ...,ip 5 ,vi = 0,F,G are of 5th degree with 
all five roots equal and G(/x) = — F(/z). 

5 

C) For ^ = Jiipi(x — Ci,y — C{), V\ =0,F,G arc of 5th degree and G(/x) = 

i=l 

—F(p), it was found that the only condition is that {cj} are roots of F. 

6 The basic cases 

6.1 Case 17: The generating system with algebraic ip 

Here we try the first term of the sequence (67): 

V> = 2 JVV (73) 

one is lead by equations (70b) and then (68) to rational expression for F 

F(A) = a 5 A 5 + a 4 A 4 + a 3 A 3 + a 2 A 2 + aiA + a + G(ji) = -F(ji) (74) 

A 
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The Lagrangian is given by 



2 .-.2 



L = J(A-m)[ ' " 



2^ F(X) F(n)' 

Jy/-F(\)FM A ■ fj_ 

(A- M )VV [ f(a) f» w 
+ 5° 6 J ^ - ^ + i (75) 



Oj, 6, c being integration constants. The second integral is 

' 2 2 

1 /xA A/t 

7 = o( A -A*)[- 



2 V _F(A) F( At ) J 

Jy/-F(A)JV) A ■ /x . 

(a- M )vv [ fw n^r 

+ lj2{ ^ + . +VMA + M) + a4]} + c A^ (?6) 

The Gaussian curvature of the configuration space carrying (75) is 

k = -i[a 5 (3A 2 + 3m 2 + 4A/x) + 2a 4 (A + fj,) + a 3 - ] (77) 
4 A xr 

In order that the system described by the Lagrangian (75) was real, it is 
necessary that the variable A varies on an interval on which F(X) > and /z on 
one of the intervals where F(fi) < 0. For the quadratic part of the kinetic energy 
to be positive definite the factor A— \i must be non-negative. The interval chosen 
for A should lie on the right of that for ii.The number of working combinations 
of intervals depends on the number of real roots of F. Figure 4 shows possible 
choices of intervals for A and \x in the case of five real distinct roots for a 5 > 0. 
When a 5 < the intervals are interchanged and some of them are excluded to 
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sustain the above conditions. 



Figure 4: Possible choices of intervals 
The possible 9 combinations are listed in the following table 



a 5 > < 



a 5 < 



A 




I 




II 




III 


in 


i 


II,III 


ii 


III 



Table III 



The general picture changes when F has equal roots or has pairs of complex 
conjugate roots. Degenerate cases will be treated in the same way. 

In a next section we shall point out an important application of the system 
I, Clebsch's case of motion of a solid in a liquid obtained as a special case. Here 
we examine integrable motions on spaces of constant Gaussian curvature. In 
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fact, when b = 0,5 = 04 = the Gaussian curvature (77) is constant n = — 
and the Lagrangian takes the form 

/ = J(A-M)[— — - - 1 



2 a 3 X + a 2 X +ai\ + a a 3 fi 3 + a 2 [i 2 + a x ^i + a Q 

y/ a 3 \ 3 + a 2 A 2 + ai A + a ^-(a 3 [i 3 + a 2 ^ 2 + ai^i + a ) 



-J 



(A - 

AA /i/i 



«3A + a 2 A 2 + aiA + a «3M + a 2A* + + a o 
+£(+*) (78) 

For this system one may use the substitutions A = Pi(£),,u = Pi (277), where Pi 
is the Weierstrass elliptic function generated by the inversion of the integral 

[ x dX 
J y a 3 X' +a 2 A +aiA + a 

to transform to isometric coordinates. However, this procedure is not transpar- 
ent enough to isolate cases of real Lagrangian systems. Suitable substitutions 
vary according to the number of real roots of the cubic polynomial, their rela- 
tive order and to the intervals on which A and fi vary. We consider somewhat 
further the case of three different real roots a,b,c (a > b > c), say. Noting the 
symmetry of the structure of the Lagrangian with respect to the variables A, ^, 
and to guarantee that the kinetic energy is positive definite, there is no loss in 
generality in assuming that A — \i > 0. The variables A, fi must also take their 
values on intervals, on which the cubic function has different signs. 

1 . When a 3 > we have three working combinations of intervals for which the 
Lagrangian describes a motion on the pseudo sphere (or on the hyperbolic 
plane H 2 ): 

(a) i. A e [a, 00), yu e [b, a] 
ii. A e [a, 00), /i e (— 00, c] 
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iii. A e [c, 6], /U e (-00, c] 

2. When 03 < we have only one working combination A G [6, a],fj, € [c, 6]. 
This leads to a motion on the sphere, which we will discuss in more detail 
in the next subsection. 

6.1.1 An integrable motion on a sphere 

In the case 2 of positive constant curvature, without loss of generality one can 
set 03 = —1. The Lagrangian and the integral can then be put in the form 



• 2 



1 r\ u A , A 



2 



2 V ^ n {a- X){X-b){X-c) (a - - /i)(/x - c) J 
J V(a - A)(A - b)(A - c y( a - - - ~c) 
4 (A-/i)VV 

AA /i/i 



x 



and 



L (a - A)(A - 6)(A - c) (a - n){b - m)(m - c) J 



7 = ( A -A*)[t TT7T TT7T 7 + 



(a- A)(A-6)(A-c) (a - fj,)(b - m)(m - c) J 

J 



- A)(A - b)(X - c)y/(a - /x)(6 - m)(m - c) x 
AA jiji 



L (a - A)(A - b)(X - c) (a - fi)(b - - c) J 
_h^X+j^ + aJ^ 
Xfj, Xfi 

In the last system it is quite easy to recognize that A, fi are the conical coordi- 
nates on the unit sphere. They are related to 3D Cartesian coordinates through 
the relations 



/ (q-A)(a-M) / (A-jKb-M) / (A-c)(/x-c) 

V (a-6)(a-c)' y V (a-6)(6-c)' Z \/( a - c )( 6 - c ) 



34 



Note also that 



a b c abc 



and 

{b + c)x 2 + {c + a)y 2 + (a + b)z 2 = X + ^ (84) 
so that A, fi are the solutions of the quadratic equation 

2 2 2 

C 2 - [(b + c)x 2 + (c + a)y 2 + (a + b)z 2 }( + abc[^- + V — + *—\ = (85) 
The potential term in the Lagrangian (80) 

Vo = ~ = (86) 

4 V 4a& C (^ + ^ + f) 

is a separable potential. The Lagrangian (80) thus generalizes the separable 
system by the introduction of the gyroscopic terms. Comparing the Lagrangian 
(80) to the results of our paper [18] we conclude that the first gives a new in- 
tegrable case of the motion of a particle on a sphere or, alternatively a new 
integrable case in rigid body dynamics under forces with scalar and vector po- 
tential. 

6.1.2 Systems on a flat manifold 

The condition of vanishing Gaussian curvature k — leads to the plane case 

v Cat (a -a) 

(a — a) x 2 — ay 2 + aa (a — a) 

= l;Ci(a-a)[- — - 2 -] 3 / 2 (87) 

2 L (a - a)x 2 - ay 2 + aa{a - a) 1 

discussed firstly in section 3.4 of our work [17]. A case presented in [23] can 
be easily obtained as a limiting case of this system. An equivalent system was 
discussed in [26] where one more limiting case is pointed out. 
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6.2 Case 17: tp is a linear function 

The choice ip = J(X + /i) leads to the expressions 



F(\) = a e X 6 + a 5 X 5 + a 4 A 4 + a 3 A 3 + a 2 A 2 + aiA + a , G(p) = -F(n) (88) 



1 Q u A A 2 i JV-F(A)F(^) r A A 



2 V \F(A) F(^) J A- M L F(A) F( M ) J 

+ i J 2 [a 6 (A + m) 3 + a 5 (A + fi) 2 } -K 1 {X + n) (89) 



and the integral 



1 AA 2 ^/i 2 

7 = 2 {x -^ [ F(x)~np) ] 

+ \j 2 {(X + M)[«e(A + m)(A 2 + /i 2 - A/i) + a 5 (A 2 + m 2 )] 

+a 4 (A + M ) 2 + a 3 (A + M )}-2if 1 A A i (90) 



/s = - J [2a 6 (A + /i) (2A 2 + 2/i 2 + A/x) 

+a 5 (3A 2 + 3^ 2 + 4Am) + 2a 4 (A + /i) + a 3 ] (91) 



6.2.1 The second case of Clebsch of motion of a solid in a liquid 

Under the restrictions a@ = = a 4 = 0, the Gaussian curvature becomes 
k = — In order that the configuration manifold was an ordinary sphere one 
may set a 3 = — 1 (say). The cubic polynomial must have three different real 
roots a, b, c(a > b > c) and the variables A, fi should be confined to the intervals 
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[b, a] and [c, b], respectively, so that A — [i > 0. We obtain 



L = t,(A-m)[ 



2 V ^ /L (a- A)(A-6)(A-c) (a - /i)(6 - ^ - c) J 
y/( a - A)(A - b)(X - c)y/(a - -//)(//- c) 
(A -A*) 

A A , 

+ ■ 



(a — A)(A — b) (A - c) (a - - - c) J 
+/ii(A + /x) (92) 

Comparing this to formulas of subsection (6.1.1) above, we find that the poten- 
tial is 

V = h 1 (X + l i) 

= /ii[(fe + c)a; 2 + (c + a)y 2 + (a + b)z 2 } 

= const - h x (ax 2 + by 2 + cz 2 ) (93) 

Comparing potential and gyroscopic terms to the list of integrable motions on a 
sphere in [18], one can show that we are dealing with the famous case of motion 
of a rigid body in a liquid named as the second case of Clebsch, characterized 
by a spherically symmetric inertia matrix. We shall not reproduce the integral 
in this case, which is well known. 

6.3 Case 19: ip = J(\-fi) 2 
For this family of systems 



F(X) = a 3 A 3 +a 2 A 2 +a 1 X + a , 
G(n) = -a 3 u 3 + b 2 /J 2 + bin + b o 

2 



" 2 HjJ^-jJ^ 



F(X) 



1 



J 2 [a 3 (A- Ai ) 2 + (a 2 + 6 2 )(A- M )] 



(A - p) 



(94) 



(95) 
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i 



+ ^ 2 (A - m)MA 2 - n 2 ) + a 2 \ + b 2 [i + at + b{\ + 



1 T . X + fi 



2 L \-n 



and 



K= 4 



-a 3 + 



2(a 2 + 6 2 )Aa* + (ai + &i)(A + /i) + 2(a + 6 ) 



(A — A*) 3 



(96) 



6.4 Case 20: ip = J^g 



For this case we get 



F(A) = A 2 (a 3 A 3 + a 2 A 2 + aiA + a ), G(m) = {i 2 (b 3 /j, 3 + b 2 /J. 2 + b\[i — a ) (97) 



Note that F(0) + G(0) = 0. After some tedious calculations we obtain as final 
result of (68) 



16AtI> - (a 3 + b 3 )ip 3 
16 a/i/- 2 - 4 J 2 



^(a 3 - & 3 )V> 2 



The Lagrangian of the present system has the form 



MM 



1 / F(A)G( M ) aa 

2 J V aV ^(a) ^(m) j 



^ A + M t2 (A + m) 2 / Q3 ^ 
A- M 8(A- M )> A j 

This Lagrangian admits the integral 

7 = 2 (A -^ )[ nA) + GR ] 



MM 



-2^-M)W^3^[^y-^yj 



2A/i — a(A + /1) 



X — fx 



/i A + /i A 



(+*) 



(98) 



(99) 



8 /i A A — \l /j, 



[a 3 -(A-3/i)-6 3 ^(3A-M)]} (100) 
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The system (99) is completely new. The fifth degree polynomials F, G have 
an identical double zero root but the other three are independent for each poly- 
nomial. The Gaussian curvature corresponding to (99) can be put in the form 

k = ~ [a 3 (3A 2 + 3n 2 + 4A^) + 2a 4 (A + ^) + a 3 ] 

(a 3 + & 3 )?/ 4 (5A - 3/i) + (02 + b 2 )y 3 (2X - M ) + (oi + h)y 2 (3X 

4(A- M ) a 



6.5 Cases 21 and 22: -0 = ^(v^i a/a*) 2 



As in the previous section, for the present two cases we obtain in final form 
L = i(A-M)[/ ^ 1 



2 V _F(A) F(^) J 



(VA±Vm) l -F(X)F( i i) Vaa 
(A-m) V A/i [ F(A) + F( A1 ) J 

~7fZ-F + ^V(Va ± vTi) 2 (102) 

VA± y^U 2 

F(A) = \(a 3 \ 3 + a 2 \ 2 + ai X + a ) (103) 



This system will have the integral 

• 2 2 

1 /xA A/i 

J = s(A-wh 



2 V \F(A) F(^) J 



7=^7= + ^ 2 [« 3 (A + m) + « 2 ](VA ± 

VX±^/JI 2 



(104) 



The two systems just obtained are different, in the sense that the difference 
in sign is significant. In fact, one can see that as A,/i approach one and the 
same root of F (different from 0) the potential term yr^^- rcmams finite for 
the positive sign but becomes infinite with the negative sign. The present two 
families of systems are new. Both families share the same configuration space. 
The Lagrangians can be explicitly in terms of global isometric variables £, r\ 
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using elliptic functions of complementary moduli. The Gaussian curvature of 
the configuration space can be calculated as 

K = -J[2a 3 (A + /x)+a 2 ] (105) 
When a 3 = this space becomes of constant curvature. 

7 Variations of cases 17 and 18 

It was noted in [16] that if we choose F(X) as the polynomial part of (74), i.e. 

F(X) = a 5 A 5 + a 4 A 4 + a 3 A 3 + a 2 A 2 + aiA + a (106) 
then one can take ip as a linear combination of five terms of the type (73) 

3 

provided { Vj,j = 1, ...,5} are the roots of the equation F(u) = 0. Here Jj are 
arbitrary constants which may be chosen real or complex in such a way that 
ip is real. The system constructed in this way involves 11 arbitrary parameters 
{ao, 0,5, Ji, J5}. When some of the roots {vj} are equal, some of the 
terms in (107) becomes dependent of the others and the number of significant 
parameters { Jj} is reduced by the number of extra repetition of roots. It turned 
out that the coalescence of roots can be made in such a way to preserve the total 
number of free parameters in the system. 

1. The first procedure was indicated in [16]. A root v repeated (k + 1) times 
leads to the replacement of k terms in (107) by terms of the form 

i=l 

so that we have the sequence of terms that can appear at possible coales- 
cence combinations are obtained from (67) by a shift v in both variables. 
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2. When the leading coefficient 05 of the polynomial (106) vanishes one of 
the roots goes to infinity. The same happens with every vanishing leading 
coefficient. It turned out that for every infinite root one can associate 
one of the terms in tp of the sequence obtained by expanding the func- 
tion ^/(A — — v)j v in powers of i . This gives again the sequence of 
expressions (66). Terms of this sequence appear in ip according to the 
number of infinite roots. 

In the following subsections we provide information, sufficient for clear and 
easy characterization of all the possible 19 cases with different combinations of 
coalescence of roots of the above two types. For every case we give only the 
forms of F, tp and the potential function V{tp). The Lagrangian is given by (55). 
The integral can be constructed in a way similar to that followed in the above 
cases. 

We shall identify each case by an ordered triplet expressing, respectively, the 
number of equal finite roots, the number of finite roots which are not equal to 
any others and the number of infinite roots. The first number may consist of 
two parts, when there are two sets of equal finite roots. Possible numbers are 
(2, 2) and (2, 3) only. In this classification the primitive case of five distinct roots 
will carry the triplet [0, 5, 0]. The case of five infinite roots [0, 0, 5], the same as 
F{\) = 1, leads to a case of pseudo-Riemannian metric (on which the kinetic 
energy is not positive definite) and we shall keep such cases out of the main list 
and list them separately. The 19 cases follow in the next two subsections 

7.1 Classification of the resulting cases 

With the appropriate choice of the intervals for A, ^i, the following 19 cases admit 
interpretation as motions on Riemannian or pseudo-Riemannian manifolds: 
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Among those cases one can easily recognize the three cases 33, 34 and 35 as 
cases of motion on spaces of constant curvature. The first of those accommo- 
dates spheres and pseudo spheres. In fact, if the roots Vi are real and ordered 
such that v\ > v 2 > vz, then one can choose for A and \i the intervals [f2,fi] 
and [2^3,2^2] , respectively, and put a 3 = —1 to be on a sphere. Similarly, if 
a.3 = 1,A € [^1,00) and jj, € (—00,^3] we have a pseudo-sphere. We do not 
write Lagrangians in specific coordinates adapted to each of the possible choices 
of intervals. In the other two cases only pseudo-sphere is possible. The system 
34 produces two possible Lagrangians. In case 35 the only possible Lagrangian 
can be written as 



L = 



1" 2 + " 2 ,.2 , -2x 

(u z + v j 



2 u 2 v 2 

■ 8J5 
+2 1 



+ 



J4U 2 



v 3 v(u 2 + V 2 ) 



8J5 



J 4 v 2 



u(u 2 + v 2 ) 



2\2 



+i28j; 



(u 2 + v 2 ) 



+J5 



32 J 4 - 



+ J3 



{u 2 + v 2 ) 



2Y2 



— 4J; 



(u 2 - V 2 ) 



-32 



Ji 



(108) 



2uv uv 

Here we have omitted terms in the linear part of the Lagrangian, which have 
no contribution to the gyroscopic function H 

Cases 36 and 37 are Euclidean plane cases in coordinates of elliptic and 
parabolic types. Case 36 is equivalent to the seven -parametric case given in 
section 3.1 of [17]. Case 37 can be shown to be representing in the xy— plane a 
system with 



V = Ax + By - a(x 2 +y 2 )[C + bx+ -a(5x 2 +y 2 )}, Q = 6ax + b (109) 

Trying to eliminate the parameter B by a rotation of the axes would result in 
new terms in both functions ft, V. Thus, this is a significant generalization of a 
case found in [16] by introducing the parameter B and freeing the parameter A 
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from a relation to a. The integral for this system is 

I = a [2{xy - yx) - (Aax + b) (x 2 + y 2 )] [y - C - bx - a (3x 2 + y 2 )] 

+A[y-bx-a (3x 2 + y 2 )] - B [x + y (b + 2ax)} (110) 

In the first fifteen cases of table IV the intervals on which the variables A, \i 
take their values could be chosen so that the configuration manifold become 
real Riemannian. The remaining four cases, in which the metric of the manifold 
has the single signature (1,-1) are listed as 38 to 41. We have just reduced the 
Lagrangians for those cases to real form by a substitution tp = where "J is 
real. No effort is done to reduce the cases obtained explicitly to the simplest 
form. 

8 More about interpretations 

In the previous sections we have constructed forty one different integrable sys- 
tems admitting a quadratic integral independent of Jacobi's integral. The pres- 
ence of several parameters in the structures of those systems gives a wide pos- 
sibility to adapt those parameters to fit concrete applications in mechanics and 
physics. So far, we have identified some possible cases of motion in the plane or 
developable surface, sphere or pseudo-sphere (or hyperbolic plane) by analyzing 
the Gaussian curvature of the metric of the system under consideration. The 
question arises, how to apply the above results to a given problem in mechanics. 
The first and principal step is to identify the metric of the configuration space 
of the given system with that of one of the systems constructed in sections 3 
- 7. This may require the use of certain transformation of the generalized co- 
ordinates of the original mechanical system to reach an isometric form with a 
coefficient A of one of the types (18) - (20) and try to adapt the candidate in- 
tegrable system by imposing some restrictions on its parameters to fit it to the 
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given one. Once this step has succeeded, we are already dealing with an inte- 
grable case of the given system and it only remains to characterize the potential 
and electromagnetic forces in in terms of the original coordinates of this sys- 
tem. If the original given mechanical system is obtained from one of more than 
two degrees of freedom by means of reduction in the sense of Routh, the cyclic 
constants of the motion would appear in the linear and potential terms in a 
certain specified form and the parameters of the candidate system may allow to 
accommodate those constants either completely or partially. Extra-parameters 
not coming from cyclic constant may then be considered for further explanation 
of their origins. In the rest of this section we indicate certain applications of 
the above results. For size considerations only some of them are considered in 
some detail and others are pointed out. 

8.1 Application to rigid body dynamics 

As an example, we try here to find all cases of the above systems, relevant 
to the dynamics of a rigid body moving about a fixed point. Reduction of 
the most general problem of motion under the action of potential forces and 
forces of electromagnetic type was reduced to the form involving a Liouville- 
type metric in [16]. We give the final form of the Lagrangian for the case of 
tri-axial ellipsoid of inertia. Let I =diag(A, B,C ) be the inertia matrix of the 
body in the principal axes of inertia at the fixed point, u> =(p, q, r) its angular 
velocity, 7 =(7 1; 72,73) a unit vector fixed in space and V a (7) , m (7) be the 
scalar and vector potentials of forces acting on the body. The Lagrangian of the 
system has the form 

L = • bj + m • w— Vq (HI) 



46 



The Euler-Poisson equations of motion can be written as 

dV 

dbl + ljx(vI + M) = 7X— — , 71^x7 = (112) 

Of 

where 

M=A(m. 7 )-(|-.m) 7 (113) 
This system admits the cyclic integral 

(wI + m)-7=/ (114) 

corresponding to the angle of precession tp around the vector 7 and allows 
reduction by the Routhian procedure to a system of two degrees of freedom. 

8.1.1 Reduction in the case of a tri-axial body 

To describe the position of the current point on the configuration space of the 
reduced system in this case, which is the Poisson sphere \j\ = 1 we shall use the 
same variables A, \i as in ([16]). Those variables are equivalent to elliptic coor- 
dinates on the ellipsoid of inertia with the standard metric. They are solutions 
of the quadratic equation 

DX 2 - NX + ABC = (115) 

N = AiB + Chl+BiC + A^l + CiA + B^l 

D = AjI + BjI + CjI (116) 



We get 



2 .2 



R = \abc{x-h)[ — a ^ — - + - 2 £ 



2 ' ' ^' l X 2 {A-\){X- B){X-C) V 2 (A- fi){B - -C) J 

+VABC[P—, x Q t ] 

X^(A-X)(X-B)(X-C) ^(A-MB-Mn-Cy 

-V (117) 
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where 



V=\[V + ^pl ] (118) 



We do not express P, Q here explicitly as lunctions ol A, /i and only compare 
the function 

i , pip , an 

n = ^=={^(A »)(B M )( M - O— \^{A A)(A - B)(A - C)-£] 

= ^ 4^5^ {/ [A + B + C ~ 2(A + M)] + S} (119) 

in which / is the cyclic constant and S is a function depending on vector- 
potential terms 

f) T 1 1 

(120) 



— 7 x( 7 I x m) 



8.1.2 Integrable cases of a tri-axial body 

We now look for cases of the constructed systems, which can accommodate 
the metric in (117). Note that here the functions F and G are fifth-degree 
polynomials with a double zero root satisfying G{pi) = — a condition 
common to most cases, but not all, of the previous sections. We find that: 



1. The case of 6.1 under the restrictions a\ = a = b = 0.a 2 = l,a 3 = 
-(1/A + 1/B + 1/C), a 4 = 4±i±<Z, a 5 = 

2. Case 24 of table IV under the restriction z/ 4 = 0. 

3. The case of 6.2 under the same restrictions, but with instead of b. 

4. The case of 6.4 under the restrictions ao = — bo = 1, ai = — b\ = — (1/A + 
1/B + 1/C),a 2 = -b 2 = ^±§±2, a 3 = -b 3 = 

It was shown in [16] that case 1 is a case of motion of a multi-connected 
tri-axial solid in a liquid, known after Clebsch [42] and case 2 is the case due 
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originally to Steklov [43] in its generalized form due to Rubanovsky [44] and 
involving circulations. Case 3 was also pointed out as a new case. The reader 
is referred to [16] for full detail. 

The last case turns out to be a new one different from known cases of rigid 
body motion. For this case one can find 



2ABC^N 2 - AABCD 
m = ^^(A(B + C) 7l ,B(C + A) l2 ,C(A + B) 73 ) 

M = -Jf(^,^,^) (121) 

where Equations (112) for this choice admit the conditional integral 

I = \{Ap 2 + Bq 2 + Cr 2 ) K{ Pll + q l2 + r l3 ) + -j=JL== (122) 

In fact, differentiating / with respect to t in virtue of (112) and using the 
geometric and cyclic integrals 

71 +72 +73 = 1, 
Apji+Bqjz + Cnz + ^^N = f (123) 



wc find 



^ = AkfC(A B){B -C)(C- A) lll2l l (124) 



and hence / vanishes identically in two cases: 

1. When fc = 0, and this reproduces the original Steklov case which is a 
general integrable case for all initial conditions. 

2. When k ^ 0, / = and this gives a new conditional case valid only on the 
zero level of the cyclic integral (123). 
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8.2 On time-reversible integrable problems in rigid body 
dynamics 

For the system (117) to be time-reversible one must have m = 0, / = 0. The 
Lagrangian becomes 

1 A 2 ' 2 1 

R = -ABC(\-n)l x2(A A)(A _ B){x c) + fi2{A _ c) ]— 4 Vo 

(125) 

It is clear from the consideration of section 1 that the existence of a quadratic 
integral is associated with the separation of variables. Comparing the last sys- 
tem to (22) and using (115) we note that it admits a quadratic integral if and 
only if 

Vb = D ^ N + ^ + MN-Vl3) (126) 
VP 

where (3 = N 2 — AABCD and v\,v 2 are two arbitrary functions. The explicit 
form of the integral in the Euler-Poisson variables is 

/ = 1(V + Bq 2 + Cr 2 ) + {N + ^ )VI{N + ^) + ( N - vgMiV - V?) 

2 vT? 

(127) 

This can be verified by direct calculation. This result was obtained in [30] 
and independently in equivalent form in [51]. In [49] an infinite sequence of 
potentials {U n ,n = l...oo} that admit quadratic integral, using symbolic com- 
putation was found. It can be easily shown that this sequence corresponds to 
the choices {v\(x) — —v 2 (x) = x n } . This is exactly a case discussed explicitly 
in [51]. Moreover, if in (126) we take vi(x) = —v 2 {x) as a full power series we 
obtain the potential 

00 

V = J2 A nUn (128) 
1 

admitting the quadratic integral. This could be expected from the separation 
property overlooked in [49] . It should be mentioned here that the authors of [49] 
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have also argued that the sequence they obtained is different from another one, 
obtained earlier by Wojciechovski in 1985 [50]. However, the latter sequence as 
provided in [49] is not separable for the tri-axial rigid body and is inconsistent 
with the existence of a quadratic integral for rigid body dynamics. 

8.3 The problem of motion of a particle on a smooth el- 
lipsoid 



Some mechanical systems do not have much importance in their own as much 
as in other problems related to them, for which they serve as geometrization. A 
typical example of a useful model in mechanics is the problem of motion on a 
smooth fixed surface. The motion on a sphere is used for modeling motion of a 
rigid body [18], in the study of B-phase of the superfluid 3 He, in the construction 
of certain wave solutions of the Landau-Lifshitz non- linear equation (e.g. [46]) 
and in the treatment of Dyson's fluid dynamical model of spinning gas clouds 
maintaining ellipsoidal shape [45] . Recently, the motion of a particle on a smooth 
tilted cone was found to be a generalization covering the problem of motion of 
the swinging Atwood machine as a special case [47]. 

The ellipsoid is one of the favorite models. It is closely related to the dy- 
namics of rigid bodies. Let a particle of unit mass be moving on the ellipsoid 

< i29 » 

under the action of forces with potential V and vector potential 1. The La- 
grangian for this problem is 

L = h 2 + l-r-V (130) 
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We shall only express f 2 in terms of elliptic coordinates u, v on the ellipsoid 

(131) 



_{a 2 - u){u - b 2 ){u - c 2 ) {a 2 - v)(v - b 2 )(v - c 2 ) 

Here the functions F and G are rational (cubic to linear). One can immediately 
notice that case I (6.1), in which F is rational (degree 6 to one) gives a case of 
motion on the ellipsoid if coefficients are suitably adapted. This case was found 
earlier in our work [19] by a different method. 

8.4 Geometrical interpretation of systems of the type (17) 

Using Hamilton's principle in the form of Jacobi (e.g. [34], [35]), one can de- 
rive the equations of the trajectories of motion of the system described by the 
Lagrangian (17a) from the variational problem 

8 J \l 2k{h - V){di 2 + dri 2 ) + Pdi -Qdr 1 = (132) 

This means those trajectories are geodesies of the Randers metric 

ds = ^2A(h-V)(d£ 2 + dr] 2 ) + Pd£ - Qdi] (133) 

Randers metrics were first studied by physicist, G. Randers [36] from the 
stand point of general relativity. Since then, many geometers have made efforts 
in investigating the geometric properties of Randers metrics as an important 
and rich class of Finsler metrics and to explore their prospective applications 
(see e.g. [37] - [40]). 

Every integrable mechanical system of the type (17) corresponds to an inte- 
grable geodesic flow of the Randers Metric (133) admitting the integral 

I = A(h- U)(|) 2 + y/2A(h-V)(P^ + Q^) + R (134) 



in which dl = \J d^ 2 + drf. 
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8.5 The Fokker-Plank systems 



Stochastic dynamical models described by Fokker-Planck equations, in the limit 
of weak noise, can be formally associated with Hamiltonian dynamical systems. 
The Hamiltonian consists of a kinetic energy quadratic term and terms linear 
in momenta , with zero scalar potential term (e.g. [9], [52]). The study of such 
Hamiltonian systems has received some attention from both points of view of 
integrability [53] - [56] and existence of invariant relations [57] . 

Gauge transformation of the Lagrangian was also used by Hietarinta [54] to 
construct integrable Fokker-Plank equations from known integrable Hamiltonian 
systems. The Hamiltonian function corresponding to the Lagrangian (17a), with 
a gauge term dZl f^ , has the form 

2 / ir? \ 2" 



Pi 



dz 



_ dZ 



v = o 



(135) 



where p^,Pri are the momenta conjugate to the coordinates £, r\ respectively. An 
integrable system (17) generates a Fokker-Plank system if there exists a function 
Z(€,r)), satisfying 



1 

2A 



P 



dZ 



dZY 



V = 



(136) 



The results we obtained for this problem are promising and will be presented 
elsewhere. 



9 Conclusion 

So far, we have made the most comprehensive analysis of the problem of con- 
structing 2D conservative Lagrangian systems involving terms linear in the ve- 
locities, which admit a complementary quadratic integral. The problem under 
consideration has proved unexpectedly rich. We have constructed and tabulated 
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41 major several-parameter families of intcgrablc systems of this type. In most 
cases the parametrization of the system can be chosen, usually in more than 
one way, so that the configuration manifold becomes Riemannian or pseudo- 
Ricmannian. The remaining few cases are real only on pseudo-Riemannian 
manifolds. 

Some of those families or their special cases or degenerations were listed in 
our previous works but the majority are either listed for the first time or involve 
more parameters in their structure. Usually, we give the Lagrangian and the 
second invariant and occasionally the Gaussian curvature of the configuration 
space, which we use for interpretation of simple results as cases of motion on 
a plane or a space of constant curvature. In application to the problem of 
motion of a rigid body under various circumstances, all known integrable cases 
of the type under consideration for this problem are restored and a new one 
is added. A systematic analysis of the results of this paper is still needed to 
isolate all possible subcases, which admit a preassigned physical or mechanical 
interpretation. Many of them are tentatively candidates for the study of motion 
of a charged particle under the action of spatially periodic electric and magnetic 
fields on various 2D spaces. 

It is worthy to say few words about the methods we used throughout the 
paper, in order to make clear their limitations. The difficulty in solving the 
problem stems from the fact that we deal with a linear PDE (52) for ip(A,[i) 
coupled with a differential condition involving nonlinearly the same function 
ip and linearly three functions of a single variable each F(A), G(/i), V(ip). This 
system is so involved that we may not be able to tell at present how far we are 
from its general solution. As we see from the last sections, F and G turned out 
to be rational in one case and polynomial of degree less than six in all others. We 
see from tables 4 and 5 how the degree of F(X) is interlinked with the number 
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of terms taken from the sequence of functions {?/>„> V'nl- It seems unlikely that 
F can be a polynomial of degree six or higher. However, we are working on 
few cases in which F, G are algebraic functions. Using certain uniformizing 
transformations, one can reduce F, G to polynomial form, but this changes the 
form of the metric coefficient (A — /x). This will be pursued elsewhere. 

In the more compact case considered in 5.3.1, we deal with simultaneous so- 
lutions of the linear PDE (52) and the nonlinear one (70a) in the same unknown 
function ip. It is unlikely that those equations admit other solutions than the 
list of table II. However, we have no proof yet that no other solutions exist, for 
example, as full infinite series of {ip n ,ip n }. 

Acknowledgement: I am indebted to an anonymous referee, whose re- 
marks helped me reorganize this work and to reformulate the methods and 
results in a much clearer way. 
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